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ABSTRACT 

We study the polarized infrared emission by Polycyclic Aromatic Hydrocarbons (PAHs), when 
anisotropically illuminated by UV photons. PAH molecules arc modeled as planar disks with in-plane 
and out-of-plane vibrational dipoles. As first pointed out by Leger (1988), infrared emission features 
resulting from in-plane and out-of-plane modes should have orthogonal polarization directions. We 
show analytically how the degree of polarization depends on the viewing geometry and the molecule's 
internal alignment between principal axis of inertia and angular momentum, which gets worse after 
photon absorption. Longer wavelength features, emitted after better internal alignment is recovered, 
should be more strongly polarized. The degree of polarization for uni-directional illumination (e.g., 
by a star) is larger than for diffuse illumination (e.g., by a disk galaxy), all else being equal. For PAHs 
in the Cold Neutral Medium, the predicted polarization is probably too small to distinguish from the 
contribution of linear dichroism by aligned foreground dust. The level of polarization predicted for 
PAH emission from the Orion Bar is only w 0.06% at 3.3 /im; Sellgren et al. (1988) report a much 
larger value, 0.86±0.28%, which suggests that the smallest PAHs may have moderately suprathcrmal 
rotation rates. Future observations of (or upper limits on) the degree of polarization for the Orion 
Bar or for dust above edge-on galaxies (e.g., NGC 891 or M82) may constrain the internal alignment 
of emitting PAHs, thus providing clues to their rotational dynamics. 
Subject headings: ISM: dust, extinction — ISM: general — infrared: galaxies 



1. INTRODUCTION 

The strong infrared emission features at 3.3, 6.2, 7.7. 
8.6, 11.3 and 12.7 /xm have been attributed to vibra- 
tional modes in pla nar Polycyclic Aroma t ic Hydrocar- 



bons (PAHs ) (e.g., Leger fc Puget 1984 Allamandola 



et al. 1985 ) . Additional strong f eatures at 16.4 and 



T7/zm (e.g., |Smith et al. 20071 have also been at- 
tribut ed to PAHs, although the identification is less cer- 
tain. Leger (1988) noted that a planar PAH molecule 
may emit partially polarized light if anisotropically il- 
luminated by a source of UV photons. The basic rea- 
sons are the following: i) UV absorption is favored if the 
molecule faces the illuminating source; ii) spinning of the 
molecule around its angular momentum preserves some 
memory of the illumination direction; Hi) the vibrational 
dipoles responsible for the IR emission features oscillate 
either perpendicular or parallel to the molecular plane. 
In particular, the C-H stretching mode (3.3 /im) and the 
in-plane C-H bending mode (8.6 /im) oscillate parallel to 
the grain plane, whereas the out-of-plane C-H bending 
mode (11.3 and 12.7/xm) oscillates perpendicular to the 
molecular plane. The strong emission features at 6.2 and 
7.7 firn are believed to arise from in-plane C-C stretch- 
ing and bending modes. In-plane modes (3.3, 6.2, 7.7, 
8.6 /im) and out-of-plane modes (11.3, 12.7 ^m) s hould 
exhibit orthogonal polarization angles (Leger 19881, and 
their electric field vector should be respectively perpen- 
dicular and parallel to the plane-of-sky projection of the 
illumination direction. 



pellgren et al. ( 1988 1 have searched for linear polariza- 
tion of the 3.3 and 11.3 /im emission features in a variety 



of astronomical sources where PAH emission is observed 
offset from the illuminating source. Their upper limits on 
the degree of polarization are of the order a few percent. 
At one position on the Orion Bar they measure a linear 
polarization of 0.86 ± 0.28 % in the 3.3 /jm feature, with 
the polarization angle consistent with being orthogonal 
to the line betw een the nebula and the star, as predicted 
by Leger (19881. However, as we show below, the polar- 



ization they report for the 3.3 /im feature is much larger 
than expected. 

In this work we present analytic formulae for the de- 
gree of polarization of the PAH emission features, when 
the emitting grains are anisotropically illuminated. We 
model PAH molecules as planar disks with in-plane and 
out-of-pla ne vibrationa l dipoles. We extend the calcu- 
lations by |Leger" (1988) to allow for an arbitrary degree 
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of disalignment between the molecule's principal axis of 
inertia ai (perpendicular to the molecular plane) and its 
angular momentum J. We discuss both the case of a 
point- like illuminating source, which may be applied to 
reflection nebulae like the Orion Bar, and of an extended 
source (e.g., a disk galaxy), which may be relevant for 
dust above NGC 891 or M82. The level of polarization 
is sensitive to the angle between the line of sight and the 
illumination direction, and to the degree of alignment 
between kx and J. Measurements of the degree of polar- 
ization can therefore provide insight into the rotational 
dynamics of PAHs. 

This work is organized as follows: in f|2]we describe our 
model for UV absorption and polarized IR emission by 
PAH molecules, commenting on uncertainties regarding 
the alignment of ai with J; in [J3]we present our results, 
both for a star-like illuminating source and for an ex- 
tended galactic disk. The reader interested primarily in 



2 



the observational implications of our work may wish to 
skip f|2]and ^3] and proceed directly to ^ where we sum- 
marize our findings and discuss how future polarization 
measurements may constrain the geometrical and rota- 
tional properties of PAHs. 

2. MODEL FOR POLARIZED EMISSION FROM PAHs 



2.2. Cross section for absorption of starlight 
The absorption cross section f or incident light with 



As discussed by Leger (1988), planar PAH molecules 
may emit partially polarized light as a result of 
anisotropic illumination by a source of UV photons. UV 
absorption is favored if the molecular plane is perpen- 
dicular to the illumination direction. Following UV ab- 
sorption, in-plane and out-of-plane vibrational modes are 
excited, producing the observed IR emission features. 

The grain angular momentum J stays approximately 
constant duri ng the whol e process of UV absorption and 
IR emission ( |Leger||l988| ): first, the angular momentum 
contributed by the absorbed UV photon or removed via 
vibrational IR emission or rotational radio emission is 
small compared to the mean angular momentum of in- 
terstellar PAHs; secondly, collisions of the emitting grain 
with interstellar atoms or ions hardly occur during the 
few seconds of IR emission; finally, Larmor precession 
of J around the interstellar magne tic field takes much 
longer than the IR emission burst (Rouan et al. 19921. 
With J conserved, some memory of the source direction 
is retained and the IR emission bands will be partially 
polarized. 



2.1. 



Illumination geometry 

We adopt the system of coordinates used by |Leger 
( 1988 1 to specify the illumination geometry and the ori- 
entation of the emitting molecule (Fig. IT]). The fixed 
coordinate system (x,y,z) is centered on the emitting 
grain; the polar axis z is along the illumination direc- 
tion and x is in the plane defined by z and the direction 
n from the molecule to the observer (Fig. [l]a). In this 
frame, the grain angular momentum J has spherical an- 
gles 9 and ip (Fig. [l]&). In order to describe the posi- 
tion of the molecule with respect to J, we define a frame 
(x',y',z') centered on the grain with z' along J and x' 
perpendicular to the plane (z,z'), as shown in Fig. [T1&. 
We choose the grain axes of inertia (a!,^,^) so that 
ai is the axis of largest moment of inertia (i.e., the prin- 
cipal axis); for a disk molecule, ai is perpendicular to 
the plane of the molecule, whereas &2 and &3 are in the 
plane. Their position in the frame (x', y', z') is described 
with Eulcr's angles (Fig.jTjc): /3, nutation angle, between 
z' (or J) and ai; ip, precession angle, between x' and the 
line of nodes, i.e., the intersection of the plane (x',y') 
with the molecular plane ({12,213); </>, angle of proper ro- 
tation, between the line of nodes and a.2. The classical 
motion of a rigid axisymmetric grain is a combined rota- 
tion around its symmetry axis a! and a precession of this 
axis around the angular momentum J with constant (3. 1 
To describe the polarization of the emitted radiation, we 
define a, the angle between the illumination direction z 
and the observer direction n, and the polarization direc- 
tions u and v, respectively parallel and perpendicular to 
y in the plane of the sky (Fig. [l]a). 

1 For a perfectly symmetric grain, the moments of inertia I2 and 
73 corresponding to £12 and a3 are equal and no nutation occurs. 



electric field E is proportional to (Leger 1988) 



AcxV|(j|E-d| 



(1) 



1,3 



where the summation is over molecular states |i) and \j) 
and d is the electric dipole moment operator. The it — tt* 
electronic transitions responsible for UV absorption in 
PAHs have (j|d|i) only in the molecular plane. For a 
disk molecule, rapid spinning around its principal axis ai 
results in averaging over the angle of proper rotation (j). 
Thus, the grain absorption cross section may be written, 
for incident unpolarized light, 

,4 oc (1 + cos 2 6) , (2) 

where O is the angle between the normal to the grain 
plane (i.e., the principal axis ai, for a disk molecule) and 
the direction of propagation k of the absorbed photon. 
In other words, when a planar PAH faces an unpolar- 
ized source, its UV absorption cross section is twice that 
when it is edge-on, because both components of the illu- 
minating electric field can be absorbed in the first case 
and only one in the second. 

The angle depends on the instantaneous orienta- 
tion of the grain with respect to the illumination direc- 
tion. However, since the precession period of ai around 
J is much shorter than the time between absorption and 
emission, we can average over the precession motion (and 
the precession angle ip) for the absorption and the emis- 
sion process independently. For a point-like illuminating 
source, all incoming rays have k = z (see Fig. [l]a). The 
corresponding ^-averaged absorption cross section, for 
fixed J and a fixed angle [3 between ai and J, is 



A*{6, p3) = 1 + cos 2 9 cos 2 (3 + - sin 2 8 sin 2 f3 
= l + C(0,/3) . 



(3) 

If the illuminating source is an extended disk galaxy, 
the angle also depends on the spherical angles 8' and 
ip' of the ray direction k in the fixed (x,y,z) frame. We 
assume that the emitting molecules are above the galac- 
tic center; z is then the direction from the center of the 
galactic disk to the grains. For an axisymmetric surface- 
brightness profile, the absorption coefficient for incom- 
ing rays with polar angle 8' can be averaged over the 
azimuthal illumination angle ip': 

A gal (9,[3,9') = 1 + cos 2 8' C(8 , (3) + sin 2 8' S (8, (3) , (4) 

where we have also averaged over the precession angle 
ip. In the previous expression, C(0, (3) is the same as in 
eq. ([3]) and we define 

S(9,(3) = i(l + cos 2 60sin 2 /3+ ^ sin 2 (9 cos 2 /? . (5) 

For the sake of simplicity, in the following we assume a 
uniform-brightness disk galaxy; for a generic axisymmet- 
ric brightness profile B(9'), the absorption cross section 
in eq. Q should be convolved with B(9'). As a special 
case, we observe that if 9' = for all incoming rays, i.e. 
k = z, we recover the absorption cross section in eq. |3| 
for a point- like illuminating source. 
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Fig. 1. — Geometry of the star-molecule-observer syst em, as presented by |Lcgcr (1988). In 6), the molecule (dashed disk) has been 
moved far from the center of (x,y, z) for clarity. See j]2.1|for details. 



2.3. Cross section for polarized emission 

The cross section for PAH e mission pola rized along a 
direction w is proportional to ( Leger|[l988 1 



■y: 



3,k 



dli) 



(6) 



where \j) and \k) are vibrational states of the emitting 

molecule. We shall call Fjj, and the relative emission 
cross sections for in-plane and out-of-plane modes respec- 
tively. For out-of-planc modes, the only non-vanishing 
terms in eq. (pi involve the component of d perpendicu- 
lar to the molecule plane, and 



F 1 



(ax 



(7) 



For in-plane modes, the matrix elements in eq. (pi are 
non-zero only for the components of d in the molecule 
plane; for a rapidly-spinning disk molecule, rotational 
invariance with respect to the angle of proper rotation cj> 
yields 

Fll=l-(a!.w) 2 = l-i^ . (8) 

By averaging over the precession angle ip for the emis- 
sion proces s in dependently from the absorption (as dis- 
cussed in {2.2 \, the emission cross sections for out-of- 



planc modes with polarization along u and v are 
F^(9, ip, /?, a) = ^(cos 2 <p + cos 2 9 sin 2 <p) sin 2 (3 



- cos 2 a 



sm' 
1 
2 



! /3 



if cos 

(sin 2 <p + cos 2 9 cos 2 <p) sin 2 (3 



(9) 



- cos 2 a sin 2 £> cos 2 </5 cos 2 /3+sin 2 a C(8, (3) 



- - sin 2a sin 29 cos ip(l 
4 



3 cos 2 /3) , (10) 



where C(9, (3) has been defined in eq. pj. We have made 
use of the fact that the orientation of the angular mo- 
mentum J is constant during the IR emission burst, as 
discussed at the beginning of Sj2] 

The corresponding emission cross sections for in-plane 
modes with polarization along u and v can be computed 
from their out-of-plane counterparts by using eq. 



Fjj (9,<p,(3, a) : 



F^(9,p,(3,a) 
Fj-(O,<p,0,a) 



(11) 
(12) 



By averaging over the azimuthal angle tp of the an gular 
momentum distribution, we simplify eqs. (|9]>-([ToJ) and 
obtain 



F^(8,0,a) = S(6,/3) , 

F^-(9, f3, a) = cos 2 a S(9, 0) + sin 2 a C(0, (3) 



(13) 
(14) 

and the corresponding (^-averaged cross sections for in- 
plane modes may then be derived from eq. (Jsl) , In the 
following, we assume that the angular momentum J is 
randomly oriented in space, but in principle eqs. ([£])-( 10 1 
and eqs. (11)-(12| could be convolved with any angular 
distribution for J , e.g., if the emitting molecules are par- 
tially aligned by the interstellar magnetic field. 

2.4. Model for the alignment of ai with 3 

The probability that a symmetric grain rotates with 
angular momentum J = |J| and nut ation angle (3 (be- 
tween ai and J) may be written (e.g., Lazarian & Draine 
19971 



dP(J,0) oc f(J) exp 



sin/3d/3dJ , (15) 



where /(J)dJ is the probability that the angular mo- 
mentum £ [J, J + dJ] , -E ro t is the grain rotational energy 
and T; a is the "internal alignment" temperature which 
parametrizes the degree of alignment between ai and J. 
In the classical approximation, the rotational energy of 
a symmetric oblate grain is 



2h 



1 



^-lW/3 



(16) 



where 1\ is the largest moment of inertia, corresponding 
to ai, and 1% = J3 < ii; for a planar symmetric molecule, 
1-2 = I3 = Ii/2. For the sake of simplicity, we assume 
that all the emitting molecules have the same moments 
of inertia. 

The absorption (eqs. Q-Q) and emission (eqs. (|9|- 
( 1 2 1 ) cross sections do not depend on J but only on thi 



depend 
Therefore, 



nly on the 
we may integrate 



angl e (3 between and J 
eq. (15) with respect to J, provided f{J) is known, and 
the resulting probability distribution will be used to de- 
scribe the degree of alignment between a! and J when we 
compute the PAH polarized emission in Sj3] In the simple 
case /(J) oc S(J — J) (or, in general, if /(J) is strongly 



4 



peaked at J), the probability distribution in eq. (15) re- 
duces to 



dP 7 (/3) 



/7exp(— 7 sin 2 [3) sin/3d/3 
V / 7re~Terfi( v /7) 



(17) 



where erfi(z) = (2/y/n) J* exp(i 2 ) dt is the imaginary er- 
ror function. We have defined a dimensionless "internal 
alignment" coefficient 



J 



h 



7^ P-l 

2hk B T ia \I 2 J T ia \I 2 



(18) 



where T rot = J 2 /2Iik B . When 7 — > 00, the molecule 
principal axis tends to be perfectly aligned with the an- 
gular momentum; when 7 = 0, ai is randomly oriented 
with respect to J. 

As discussed at the beginning of Sj2] J stays approx- 
imately constant between UV absorption and IR emis- 
sion, and so does T rot ; however, the molecule internal 
alignment temperature T la may substantially change if 
part of the absorbed photon energy is transferred to ro- 
tational degrees of freedom. We account for the uncer- 
tain energy exchange between vibrational and rotational 
modes by parametrizing the grain internal alignment be- 
fore UV absorption and during IR emission respectively 
with 70 (corresponding to T; a = To) and "f r < 70 (corre- 
sponding to Tj a > To). In principle, 7 r should be different 
for each IR emission feature, since it may be thought of 
as the internal alignment coefficient when most of the 
radiation in that band is emitted. 

We now discuss another plausible choice for the 
parametrization of the alignment between ki and J prior 
to UV absorption. If many collisions with hydrogen 
atoms occur in the interval between two UV absorptions, 
and there are no other torques acting, the grain will be 
driven towards "Brownian rotation" with /(J) oc J 2 . If 
there is no vibrational-rotational energy exchange (i.e., 
the PAH acts like a rigid rotator) , the internal alignment 
temperature T before UV absorption will approximately 
equal the gas kinetic temperature T gas . Integration of 
eq. (15 1 with respect to J assuming /(J) oc J 2 yields a 
probability distribution 



dP e (/3) 



e(e~ l) 1 / 2 sin/3d/3 
2(e-cos 2 /3)3/2 



(19) 



which does not depend on T gas but only on the geomet- 
rical properties of the grain via e = — J2). 

In the following, we para metr ize the disalignment be- 
tween ki and J by using eq. ( 17 1 both before UV absorp- 
tion (with 70) and during iR emission (with j r < 70). 
However, in the Appendix we present analytic formulae 
for the expected degree of polarization if the grain align- 
ment before UV absorption can be described by eq. ( 19 1, 
but eq. (17 1 still holds during IR emission. 



2.5. Estimating 70 and 7 r 



In the dust model of |Draine fc Li| ( |2007 ), the 3 - 
13 /irn emission features are produced primarily by PAH 
molecules or clusters containing between N c w 25 and 
~ 1000 carbon atoms, and the 17 /xm complex is mainly 
due to PAHs with N c w 2000. For purposes of estimating 
rotational kinetic energies and the density of vibrational 
states, we will take N c = 200 as a representative value, 



with a volume-equivalent radius a « 7.5 A. We suppose 
that all the emitting molecules are axisymmetric and pla- 
nar, with I1/I2 = 2. 

Following UV absorption, the alignment between the 
molecule principal axis ai and angular momentum J, 
and thus the internal alignment temperature T la , de- 
pends on the efficacy of internal energy exchange be- 
tween lattice vibrational modes and rotational modes. 
The Intramolecula r Vibration-R otation Energy Transfer 
(IVRET) process (Purcell 1979[ ), due to imperfect elas- 
ticity of the molecule when stressed by centrifugal and 
Coriolis forces, allows energy exchange bet ween rotation 



and vibrations on a timescale 



10- 



s (Rouan et al. 
19921, much shorter than the duration ~ 1 — 10 s of the 



IR emission burst. This means that, while the molecule is 
cooling after UV absorption, its internal alignment tem- 
perature T a tends to be equal to the instantaneous vi- 
brational (lattice) temperature T v ;b- Following photon 
absorption, the lattice may be heated up to a tempera- 
ture T v ;b ~ 300 — 1500 K, depending on the grain size 
and the photon energy. The grain cools as IR energy is 
radiated, and we estimate T v ;b « 800 K when most of 
the 3.3 /jm emission takes place, T vib m 300 K for the 
7.7 //m emission, T v jb ~ 200 K for the 11.3 /im emission, 
and T v i b ~ 120 K for the 17 /Ltm emission, unless the lo- 
cal radiation field is so intense to prevent the lattice from 
cooling down to such temperatures. 

The internal temperature T a follows T V jb while the 
grain is cooling, with T; a « T v i b as long as the vibra- 
tional energy levels are sufficiently closely spaced to al- 
low energy transfer between vibrations and rotation. 2 
When the separation AE of vibrational levels exceeds 
~ hid Iot , the IVRET process ceases to operate and the 
rotational modes decouple from the lattice. The density 
of states can be calculated using a model norm al mode 
spectr um and the Beyer-Swinehart algorithm ( |Draine| 



& Li 



|2001| ): for N c = 200, the density of states at 
vibrational energy E/hc « 250 cm -1 is hcdN/dE w 
1/(0.18 cm -1 ). This should still allow energy exchange 
in quanta AE/hc ~ (w ro t/27rc) ~ 0.33 cm -1 for a grain 
spinning at w rot /27r w 10 GHz, so that the rotational and 
vibrational modes will decouple only when the lattice vi- 
brational energy has dropped below E/hc ~ 250cm . 
For our model PAH, this is the average energy if it were 
in contact with a ~ 65 K heat bath, which suggests 
T) ~ 65 K to describe alignment of ai with J prior to 
photon absorption. As we discuss below, this is actu- 
ally a lower limit for To, which is presumably attained in 
cold interstellar clouds, whereas in bright photodissoci- 
ation regions like the Orion Bar the radiation field is so 
dense that PAH grains stay always hotter. 

As anticipated in §2.4[ we assume that the grain dis- 
alignment between principal axis ai and angular momen- 
tum J may be described by the probability distribution in 
eq. ( |17[ ) with "internal alignment" coefficient 70 just be- 
fore UV absorption, and 7 r during IR emission. The ratio 
T ro t /Tj a depends on environmental conditions, hence so 

2 If ai and J are not parallel and I2 j£ I3, the molecule nu- 
tates. Because of nutation, the rotational motion around ai is only 
quasiperiodic, and the centrifugal and Coriolis stresses responsible 
for the IVRET process will have Fourier components over a range 
of frequencies around the mean rotation rate LJ TO t (i.e., the rota- 
tion rate averaged over the nutation period). This facilitates the 
coupling of rotation and vibrations via the IVRET process. 



5 



TABLE 1 

Selected Cases for a Disk Molecule with Nc = 200 



Parameter 


case (a) 


case (b) 


case (c) 


case (d) 




CNM 


Orion Bar 


STR 


ai || J 


Tot (K) 


80 


220 


1000 


oo 


To (K) 


65 


150 


60 




7o 


1.2 


1.5 


17 


oo 


7 r (3.3/um) 


0.10 


0.28 


1.3 


oo 


7r(7.7^m) 


0.27 


0.73 


3.3 


oo 


77.(11.3/im) 


0.40 


1.1 


5.0 


oo 


7r(17/Ltm) 


0.67 


1.5 


8.3 


oo 


pl(7r/2), 3.3/im (%) 


0.02 


0.06 


1.29 


7.69 


pI(tt/2), 7.7/um (%) 


0.05 


0.17 


3.29 


7.69 


p|(7r/2), 11.3Atm (%) 


-0.14 


-0.53 


-8.56 


-14.29 


p^(7r/2), 17Atm (%) 


-0.25 


-0.73 


-10.54 


-14.29 



does 7 in eq. (18 1. We now discuss four possible cases, 
summarized in Table [T] for each case, we compute the 
value of 70 prior to absorption of starlight photons, and 
the values of j r during emission of radiation in the 3.3, 
7.7, 11.3 and 17 /jm bands. 

2.5.1. Case (a): CNM 

In H I clouds with gas kinetic temperature T gas ps 
100 K (the "Cold Neutral M edium", or CNM), the ro- 
tational excitation model by |Draine fc Lazar ian (1998) 
predicts that a PAH with N c = 200 should be spinning 
with rotational temperature T rot ~ 80K, which corre- 
sponds to a rotation frequency w rot /27r ps 10 GHz. Since 
the molecule's angular momentum stays approximately 
constant between UV absorption and IR emission, we 
may assume T rot » 80 K during the whole process. If 
T ps 65 K prior to UV absorption, as discussed above, 
the internal alignment coefficient of our model PAH will 
be 70 ~ 1.2; the values of 7 r for different IR emission 
bands are listed in Table CD 

We note that, if gas collisions were frequent enough to 
drive the emitting grains towards "Brownian rotation" 
between two subsequent UV absorptions, then the prob- 
ability distribution in eq. (19 1 should be used instead of 
eq. ( |l7| ) to describe the internal alignment prior to UV 
absorption, whereas eq. (17 1 may still be used during IR 
emission with an appropriate choice for j r . However, as 
discussed in the Appendix, identical polarization results 
may be obtained in the current formalism if we employ an 
"effective" internal alignment coefficient 70 ~ 1.0 before 
UV absorption, similar to the value 70 = 1.2 adopted for 
case (a). 

2.5.2. Case (b): The Orion Bar PDR 

The Orion Bar, illuminated by the Trapezium stars, is 
an example of a bright photodissoci ation region (PDR) , 
with strong PAH emission at 3.3 /im ( |Tielensetal.|1993| . 
Physical condition s in the photodis sociation zone have 
been discussed by Allers et al. ( 2005 ), who estimate n H ps 
7 x 10 4 cm" 3 , T ga l « 1000 K, andx- 3 x 10 4 , where n H 
is the hydrogen number density and x is t he specific radi- 
ation energy density at 1000 A relative to the value in the 
local interstellar radiation field. For these conditions, the 



rotational excitation model of Drainc & LazaxianJ ( 1998 1 
yields a rotational temperature 1 ' YO t « 220 K for a PAH 
with N c = 200. 

In the intense radiation field of a bright PDR, the 
vibrational energy of our model PAH is always large 



enough (and the vibrational levels sufficiently closely 
spaced) that the internal alignment temperature T; a 
roughly equals the instantaneous vibrational tempera- 
ture T v ib at any time. For a PAH molecule with N c = 200 
in a radiation field with \ ~ 3x 10 4 , the lattic e can hardly 
cool below T vib ps 150 K ( |Draine fc Ll||2001J ), so that in 
a PDR like the Orion Bar 'Jo ~ 150 K may be appropri- 
ate prior to UV absorption. The corresponding internal 
alignment coefficient will be 70 ps 1.5 for T rot ps 220 K. 
Values of j r for different emission bands are given in 
Table [l] since T V jb never drops below ps 150 K, we take 
Tj a pa 150 K as the characteristic internal alignment tem- 
perature for emission in the 17 /zm band. 

2.5.3. Case (c); Suprathermal Rotation (STR) 

As discussed in §2.5. 1| and |2.5.2| PAH molecules in 
the CNM should have moderately sub-thermal rotation 
rates (T rot < T gas ), whereas T rot < T gas in bright PDRs 
like the Orion Bar. Observations of microwav e emission 
from spinning grains in the diffuse I SM (e.g., LDoblerfc 
Finkbeiner||2008| | Dobler et ai"1|2008[ ) and in PDRs (e.g., 
JyJ) appear to be consistent with these 



Casassus et al. 



estimates. 

However, it is of interest to explore the possibility that 
some PAHs might be subject to systematic torques which 
could spin them up to higher rotation rates. For example, 
inelastic collisions between a grain and hydrogen atoms 
may be followed by ejection of H2 molecules when a C- 
H bond is broken through photo- or thermo-dissociation 
triggered by the absorption of UV photons. If recombi- 
nation of H atoms preferentially occurs at a few catalytic 
centers spread over the grain surface, and the H 2 ejec- 
tion is systematically asymmetric, the resulting torque 
may significantly spin up the PAH molecule. 

For case (c) we assume T rot ~ 10T gas ps 1000 K for 
PAHs in H I clouds, corresponding to a rotation fre- 
quency w ro t/27r ps 30 GHz for our model grain with 
N c = 200 carbon atoms. For PAHs spinning at 30 GHz, 
efficient energy transfer between vibrations and rotation 
via the IVRET process is suppressed, due to insufficient 
density of vibrational states, when the lattice cools be- 
low T V ib ~ 60 K, hence we take To ps 60 K as the internal 
alignment temperature prior to UV absorption, giving an 
internal alignment coefficient 70 ~ 17. Values of j r for 
different PAH emission bands are listed in Table [l] 



2.5.4. 



ai 



Case (d); 

The physics of internal relaxation in cold, spinning 
grains is not well understood, and perhaps there is a slow 
process that couples the rotational degrees of freedom 
to the lattice down to very low temperatures, regardless 
of the lower limit imposed by the density of vibrational 
states. Since in H I clouds the interval between two sub- 
sequent photon absorptions is long (~ 10 7 s), there is a 
lot of time for a weak process to act. Such a process 
would lead to 70 3> 1, and we discuss 70 = 00 as a lim- 
iting case. 

Also, if the coupling between rotation and vibrations 
after UV absorption took place on timescales much 
longer than the duration of the IR emission burst (~ 
1 — 10 s), the internal alignment temperature during IR 
emission would stay roughly equal to the value before 
UV absorption. If 70 3> 1, this would imply j r ^> 1 as 
well. We therefore consider 70 = 00 and 7 r = 00 as the 



G 



limiting case in which ai is perfectly aligned with J both 
before UV absorption and during IR emission. 

3. POLARIZATION FROM PAHs ILLUMINATED BY A 
STAR OR A DISK GALAXY 

We now compute the degree of polarization expected 
for in-plane and out-of-plane vibrational modes of PAH 
molecules when anisotropically illuminated by a star or 
a disk galaxy. We assume, for the sake of simplicity, that 
the angular momentum of emitting grains is randomly 
oriented in space, although it may be partially aligned by 
the interstellar magnetic field. Also, we suppose that the 
illuminating disk galaxy has a uniform surface-brightness 
profile. However, the formalism developed in ^J2] can be 
applied without such restrictions, and we discuss below 
how to relax these two assumptions. The grain alignment 
between ai a nd J is described via the probability distri- 
bution in eq. (17 1 both before UV absorption (dP 7o ) and 
during IR emission (dP 7r ), with "internal alignment" co- 
efficients 70 and j r respectively. 

For a population of PAH molecules illuminated by a 
point source, the emission polarized along a direction w 
(= either u or v - see Fig. Ilia) for in-plane (||) and out- 
of-plane (_L) modes is 

(a, 70, 7r)oc/ dcos6> dip dP 7o (/3 ) / dP 7r (/? r ) 
J-i Jo Jo Jo 

(20) 



xA ir (9,(3 )F^ ± (9,ip,f3 r ,a) , 



where A+ and Fw are the absorption and emission cross 
sections computed in §2.2| and |2.3| respectively. If the 
grain angular momentum is not isotropically oriented in 
space, it is easy to incorporate the corresponding (9, ip)- 
distribution function into the above integral. The polar- 
ized emission from PAHs above the center of a uniform- 
brightness disk galaxy is 

/ fat^(a J 7o,7r,w)cx/ d cos 6 I dip dP 7Q (/3 ) / dP 7r ((3 r ) 



Mcoefl' 

■\ar.i ^2 



A^(9,Po,9')Fy(e,^p r ,a) ,(21) 



where f2 = 2n(l — cosw) is the solid angle subtended 
by the galactic disk as seen from the emitting molecules. 
A generic axisymmetric surface-brightness profile B(9') 
can be easily included in the previous expression. 

The resulting degree of polarization for in-plane and 
out-of-plane modes is respectively 

j||,J- _ r||,-L 

^-tWt^ ' (22) 

and will therefore be positive for emission polarized along 
u and negative for polarization along v. 

3.1. Polarization from PAHs illuminated by a star 



By integrating eq. ( 20 1 , we obtain the degree of polar- 
ization for in-plane and out-of-plane modes from a pop- 
ulation of PAH molecules illuminated by a star (or a 
generic point source) 



p!(a,7o,7r) : 



pf(a>7o,7r) : 



3 sin a 



640 h{lo)Hlr) + 3 cos 2 a - 1 



-6 sin a 



320 h{lo)Hlr) - 3 cos 2 a + 1 



(23) 
(24) 



where a is the angle between the line of sight and the 
illumination direction (see Fig. [l]a), and 



.91(7) = V7T7erfi(V7) , 
52(7) = 6\/7e 7 - (3 + 27) v / 7rerfi(^/7) 
Hi) = Si (7)/S2(7) ■ 



(25) 



The function ^(7) monotonically decreases from h(0) — > 
00 to h(j — > 00) = 1/4. Since h(0) — > 00, if the grain 
principal axis of inertia is randomly oriented with respect 
to the angular momentum either before or after UV ab- 
sorption, the emitted radiation is completely unpolar- 

ized. Since ^1(7) > 1/4, we have pi > and p^ < for 
any cho ice of 70, 7 r and a. According to the definition in 
eq. ( 22 1 , this means that in-plane and out-of-plane modes 
are polarized respectively along u and along v, confirm- 
ing Leger's (1988) prediction that the polarization direc- 
tion for in-plane and out-of-plane modes is respectively 
orthogonal and parallel to the plane-of-sky projection of 
the illumination direction (see Fig. 1 a). 

la d 



For fixed 70 and j r , eqs. (23 1 aa 



show that the 

IR emission bands should be maximally polarized when 
a = 7r/2, i.e., the line of sight to the emitting PAHs is 
orthogonal to the star-molecule direction. For this op- 
timal viewing geometry, Table [T] reports the degree of 
polarization expected for the 3.3, 7.7, 11.3 and 17 /im 
emission features in different environmental conditions, 
as discussed in |2.5| we have assumed the 1 7 /im band to 
arise from out-of-plane modes. Figure|2]shows the depen- 
dence of p\ (tt/2) and p+(ir/2) on the internal alignment 
coefficient 70, with the different curves corresponding to 
different values of the ratio "f r /lo- Because 7 r < 70, 
the region to the left of the solid lines is excluded. In 
most cases, the depend enc e on the viewing angle a in 
the denominator of eqs. (23 1- (24 1 can be neglected, hence 

pi! (a, 70, 7r) ~ P*(7r/2, 70, 7 r ) sin 2 a, and similarly for p+ 
(compare solid and dotted lines in Fig. |3|. 

For the extreme case of perfect alignment between ai 
and J during both UV absorption (70 = 00) and IR 
emission (7,. = 0), corresponding to case (d) in ^2.5 



eqs. (23) and (24 1 may be simplified: 



pi! (a) 



sin 2 a 



13 



cos z a 

■ 2 

— sin a 

7 — cos 2 a 



(26) 
(27) 



3.2. Polarization from PAHs illuminated by a 
uniform-brightness disk galaxy 

The degree of polarization for in-plane and out-of-plane 
modes from a population of PAH molecules above the 
center of a u niform-brightness disk galaxy is obtained 
from eq. (21 1: 



Pgal( a >70,7r,w): 



3 sin a 



1280 



(cos 2 w+cos oj) 



3 cos 2 a — 1 



pL\(a,lo,lr,u) = 



-3 sin a 



640 



(cos 2 OJ+COS uj) 



3 cos 2 a + 1 



■ (28) 
(29) 



where f2 = 27r(l — cosw) is the angle subtended by the 
galaxy as seen from the emitting grains. The limit O — > 
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Fig. 2. — For PAHs illuminated by a star, dependence of the 
polarization p + on 70, for optimal viewing geometry (a = tt/2). 
The different lines correspond to different values of the ratio 
7r/7o: 7r = 70 (solid), 7 r = O.570 (dashed), y r = O.2570 
(dotted), and 7 r = 0.1 70 (dot-dashed). The polarization values 
reported in Table [I] are shown as solid triangles (3.3 (im), open 
triangles (7.7 /^m), solid circles (11.3 /im), and open circles (17 fim, 
assumed to arise from out-of-plane modes), with different values 
of 70 corresponding to case (a) (70 = 1.2), case (b) (70 = 1.5) or 
case (c) (70 = 17). 

recov ers the c ase of a point-like illuminating source in 
eqs. (23) and (24 1, whereas for an infinite disk (O — » 
2tt) the degree of polarization tends to zero. Figure [1] 
shows, for the most favorable viewing geometry (a = 
7r/2), the dependence of the degree of polarization on 70 
for PAHs illuminated by a disk galaxy with = w; the 
different curves correspond to different values of "fo/ln 
as explained in {3.1 For diffuse illumination with f2 = it, 
we find levels of polarization that are ~ 40% of the values 
found for uni-directional illumination. 

The extreme case of grains with perfect internal align- 
ment (ai || J du ring both UV absorption and IR emis- 
sion, case (d) in £ 2.5 1 yields 



II 3 sin a 

PgaiK w ) = gOfWw + cosw)- 1 + 3 cos 2 a - 1 



—3 sin a 
40(cos 2 u) + cosw) -1 — 3 cos 2 a + 1 

4. DISCUSSION AND SUMMARY 



,(30) 
.(31) 



We have obtained analytic formulae for the degree of 
polarization expected for the 3.3, 6.2, 7.7, 8.6, 11.3, 12.7, 
16.4, and 17 /im emission features when the emitting 
PAHs are anisotropically illuminated by a source of UV 
photons. We model PAH grains as planar molecules with 
in-plane and out-of-plane vibrational dipoles. Vibra- 
tional modes oscillating in the molecular plane (respon- 
sible for the 3.3, 6.2, 7.7, and 8.6 /im emission bands) are 
polarized perpendicular to the plane-of-sky projection 
of the illumination direction, whereas for out-of-plane 



Fig. 3. — For PAHs illuminated by a star, dependence of the 
polarization £>*(«, 70, 7r) on the viewing angle a. Da shed line: 
polarization for (70, 7r) = (1 0, 2), similar t o case (c) in j ]2,5| The 
polarization values quoted by |Leger| Jl988l for a disk molecule are 
shown as solid squares, and are seen to match the values calculated 
for (70, 7r) = (10,2). Solid line: polariz ation for perfect internal 
alignment (70 , 7r ) = (00,00) (see eqs. (J26|-([27l), corresponding 

to case (d) in [2.5 The dotted line is p\' (vr/2) sin 2 a for 
(70i7r) = (00 , 00 ) , which closely approximates the dependence on 
-d27b. 



a of eqs. 
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Fig. 4. — For PAHs illuminated by a disk galaxy with Q = tt, 
dependence of the polarization p ga i on 70, for optimal viewing 
geometry (a = tt/2). The different lines correspond to different 
values of the ratio 7 r /70i as explained in the caption of Fig. [2] 
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modes (responsible for the 11.3 and 12.7 /im features) 
the polarization is along th e sour c e-mol ecule direction, 
as originally pointed out by Leger ( 1988 ) . The fact that 
the predicted polarization directions for in-plane and out- 
of-plane modes are orthogonal provides a check on the 
contribution to the polarization from linear dichroism by 
aligned foreground dust, as this would be expected to 
produce polarization in a single plane. It also means 
that polarization measurements could be used to diag- 
nose the character (in-plane or out-of-plane) of the modes 
responsible for the poorly-characterized 16.4 and 17 /im 
emission features. 

Our analytic formulae show explicitly how the degree 
of polarization depends on the viewing geometry (i.e., 
the angle a between the line of sight and the illumina- 
tion direction) and the "internal alignment" coefficients 
7o, describing the alignment between the grain principal 
axis ai and angular momentum J just before UV absorp- 
tion, and 7 r , characterizing the internal alignment during 
the brief period of IR emission following UV excitation. 
The degree of polarization is maximal for a = tt/2 and, 
at fixed a, it increases with 70 and j r , as the princi- 
pal axis and angular momentum become more closely 
aligned. We have discussed both the case of a point- 
like illuminating source, which may be applied to re- 
flection nebulae and PDRs such as the Orion Bar, and 
of an extended source (e.g., dust above a disk galaxy 
like NGC 891 or M82). The degree of polarization for 
uni-directional illumination is higher than for diffuse il- 
lumination, all else being equal; however, if the galaxy's 
surface-brightness profile is strongly peaked at the galac- 
tic center, the galaxy would resemble a point source, 
with polarization levels comparable to the case of uni- 
directional illumination. 

The value of j r is expected to increase with increasing 
wavelength of the emission feature, since such a band will 
be mostly emitted when the grain is cooler and a x and J 
are expected to be more closely aligned. Therefore, when 
comparing two IR emission features arising from modes 
of the same character (either in-plane or out-of-plane), 
the band with longer wavelength should be more strongly 
polarized. Table [l] reports, in case of uni-directional il- 
lumination and optimal viewing geometry (a = tt/2), 
the predicted degree of polarization for the 3.3, 7.7, 11.3, 
and 17 /im emission features, for four exemplary choices 
of environmental conditions. 3 

For dust in the CNM, significant disalignment between 
ai and J (70 « 1, 7 r < 0.5 70) is expected to strongly 
suppress the intrinsic polarization of the IR bands, with 
the 3.3 /im emission feature predicted to be polarized by 
only as 0.02% (for our estimated 70 = 1.2 and j r = 0.10), 
presumably too small to distinguish from polarization 
due to foreground linear dichroism. Longer wavelength 
features are expected to be more strongly polarized (e.g., 
-0.14% for the 11.3 /im feature, with j r = 0.40), but the 
polarization levels are still observationally challenging. 
These conditions may apply to PAHs above the disk of 
an edge-on galaxy, such as NGC 891 or M82; however, 
the polarization degrees quoted above are computed for 

3 The polarization results in Table [l] assume that the PAH 
angular momenta are randomly oriented in space; if the grains are 
partially oriented, e.g., by the interstellar magnetic field, we would 
expect a higher degree of polarization. 



the case of uni-directional illumination, which is realized 
only if the galaxy's luminosity is strongly concentrated 
at its center. If the galaxy appears to the emitting grains 
as an extended source, the PAH polarization will be even 
smaller. 

The Orion Bar PDR is expected to be somewhat 
more favorable for the PAH internal alignment because 
of higher molecular rotation rates, but even there the 
polarization of the 3.3 /im feature is predicted to be 
only sa 0.06% for our estimated (70, 7r) = (1 -5, 0.28). 
If the polari zation 0.86 ± 0.28% measured by |Sellgren 



et al. 



1988) for the 3.3 /tin emission feature is correct, 



it imp its that we underestimated the internal align- 
ment of emitting PAHs both before UV absorption (70) 
and during IR emission (7,.). A polarization of 0.86% 
for optimal viewing angle a = tt/2 would require, e.g., 
(70, 7r ) = (2.3,2.3), (3.0,1.8), (00, 0.81), etc.: all solu- 
tions must have j r < 70, hence 70 > 2.3 and 7 r > 0.81. 
Such good internal alignment is not expected unless the 
smallest PAHs are rotating suprathermally. The polar- 
ization at 11.3 /im is predicted to be considerably larger 
than the polarization at 3.3 /im (see Tablell]). Additional 
polarization measurements of the Orion Bar at 3.3 and 
11.3 /im a re needed. If the larg e intrinsic polarization de- 



tected by Sellgren et al.Ml988b at 3.3 /im is confirmed, it 



will call into question current thinking with regard to the 
rotational dynamics of PAHs in PDRs, since it seems to 
require at least moderately suprathermal rotation rates. 

Comparison of our analytic results with polarization 
measurements may provide useful constraints on the ge- 
ometrical and dynamical properties of PAH molecules, 
concerning in particular: i) their planarity, which has 
been assumed a priori in this work, although it is still 
uncertain, especially for larger PAHs; ii) the efficacy of 
internal relaxation processes that can exchange energy 
between vibrational and rotational modes; Hi) the possi- 
bility of systematic torques that may spin the grains up 
to suprathermal rotation; and iv) a determination, via 
the polarization direction, of the character of the vibra- 
tional modes contributing to a given IR feature, which 
will be of particular interest for the poorly-characterized 
16.4 and 17 /im bands. 

If the polarization is too small to be detected, the inter- 
pretation will not be clear: it could be the result of PAH 
non-planarity, or a mixture of in-plane and out-of-plane 
contributions to the mode (this may be an issue for the 
16.4 and 17 /im features), or (more likely) it could simply 
derive from poor alignment of the grain principal axis ai 
and angular momentum J. However, if a high degree of 
polarization is observed, it will imply both that PAHs 
are planar and that their principal axis is well aligned 
with the angular momentum. This would be valuable 
information concerning the nature of PAHs and their ro- 
tational dynamics. 

BTD thanks Chris Packham and Charles Telesco for 
helpful discussions on the possibility of polarization of 
the PAH features. We are grateful to R.H. Lupton for the 
availability of the SM graphics program. This research 
was supported in part by NSF grant AST-0406883. 
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APPENDIX 

POLARIZATION FROM PAHS IN BROWNIAN ROTATION 

We have assumed above that the internal alignment between the principal axis ai and angular momentum J of 
emitting PAHs can be parametrized via the probability distribution in eq. (17 1 both before UV absorption (with 
7o) and during IR emission (with 7,. < 7o)- However, if the interval between two subsequent photon absorptions is 
longer than the time required for the PAH molecule to collide with its own mass of gas, and if the IVRET process is 
ineffective at exchanging energy between the grain rotation and the (generally colder) vibrational modes, the molecule 
will be driven towards "Brownian rotation" with inte rnal alignment temperature To ss T g as p rior to UV absorption. 
In this case, the probability distribution dP € in eq. ( 19 1 should be used instead of eq. ( 17 1 to describe the grain 
internal alignment before UV absorption, whereas dP lr in eq. ( 17 1 with internal alignment coefficient j r may still 
be appropriate during IR emission. If the molecule angular momentum is randomly oriented in space, the polarized 
emission along a direction w (= either u or v - see Fig. Tlla) for in-plane and out-of-plane modes from a population of 
PAHs illuminated by a point-like source is 



dcos 9 



dtp dP e (p Q ) dP lr (p r )A^0,p o )Ft ± (0, l p,P r ,a) 



The corresponding expression for PAH molecules above the center of a uniform-brightness disk galaxy is 



dcos 9 



dtp / dP e (/? ) / dP lr {P r ) 



dcosf 

n 



where ft = 2tt(1 — cos a;) is the solid angle subtended by the galactic disk as seen from the emitting molecules 
degree of polarization can then be computed from eq. ( 22 1 . Let us define 



q(e) =3e[l - v , e-larcsin(l/v / e)] - 1 



(Al) 

(A2) 
The 

(A3) 



where e = — ^2)- Making use of h("f) as defined in eq. (25 1, the polarization for in-plane and out-of-plane modes 

in case of uni-directional illumination is 



pl(a,e,7 r ) = 



3 sin a 
320h(j r )/q(e) +3 cos 2 a 
—3 sin 2 a 



lfflh(i r )/q(e) -3 cos 2 a + 1 
whereas for diffuse illumination by a uniform-brightness disk galaxy we have 



pL( a > e >7r,^) = 



640 



3 sin a 

hM/q(*) I 3 cog2 ~ _ 1 



(cos 2 CJ+COS uj) 



Pgal(">e>7r,^) : 



-3 sin a 



320 



Kl r )/q(e) 
(cos 2 w+cos cj) 



3 cos 2 a + 1 



(A4) 
(A5) 

(A6) 
(A7) 



Comparison of these formulae with cqs. ( 23 1-( 24 ) and eqs. (28 l-( p9| sh ows that the polarization results obtained here 
can be exactly reproduced by the probability distribution in eq. ( |17[ ) if we adopt an "effective" internal alignment 
coefficient 70 before UV absorption satisfying 2/1(70) = l/q(e). For e = 2, as appropriate for our model PAH with 
I1/I2 = 2, the solution is 70 ~ 1.0. 



